Vortices near surfaces of Bose-Einstein condensates 
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The theory of vortex motion in a dilute superfluid of inhomogeneous density demands a boundary 
layer approach, in which different approximation schemes are employed close to and far from the 
vortex, and their results matched smoothly together. The most difficult part of this procedure is the 
hydrodynamic problem of the velocity field many healing lengths away from the vortex core. This 
paper derives and exploits an exact solution of this problem in the two-dimensional case of a linear 
trapping potential, which is an idealization of the surface region of a large condensate. It thereby 
shows that vortices in inhomogeneous clouds are effectively 'dressed' by a non-trivial distortion 
of their flow fields; that image vortices are not relevant to Thomas-Fermi surfaces; and that for 
condensates large compared to their surface depths, the energetic barrier to vortex penetration 
disappears at the Landau critical velocity for surface modes. 



INTRODUCTION 



Although much has long been known about quan- 
tized vortices in superfluids of constant homogeneous 
density, vortices in trapped dilute Bose-Einstein conden- 
sates move and interact within significantly inhomoge- 
neous background clouds. The motion of vortices in 
an inhomogeneous superfluid, including their stability 
and location in a rotating condensate, has been the sub- 
ject of several recent theoretical works, both numerical 
1 1 1111 and analytic || | |fij 0, || g| & 
After years of effort, vortices in dilute condensates 
are now also a subject of active experimental study 
II H [l7| |l| |l| H H |§. Quantitative compar- 
ison between theory and experiment therefore requires 
refinement of theoretical methods, in order to go beyond 
qualitative or order-of-magnitude estimates, and derive 
more rigorous predictions. 



This paper will take a step in this direction, by 
providing an exact solution (in a non-trivial but two- 
dimensional case) to a part of the problem that has usu- 
ally been treated with uncontrolled approximations. Our 
specific conclusions will include the verdicts that image 
vortices are inapplicable to surface effects in experimental 
condensates, and that the often-cited energetic barrier to 
vortex penetration || disappears at velocities no greater 

well be- 



than the surface mode critical velocity 23, 
low previous estimates jl[ [l3|. This conclusion applies 
in the experimentally relevant limit of a condensate large 
compared to its surface depth, and the lowering of crit- 
ical rotation frequencies which it implies is an addition 
to any effects of asymmetric potentials [O. The case 
we examine also provides a general warning against ne- 
glecting the nontrivial hydrodynamics of inhomogeneous 
superfluids. 



Boundary layer theory 

In a dilute superfluid whose macroscopic wavefunction 
^ is governed by the Gross-Pitaevskii equation 

ihd t ^ = -^H + v{x)^ + ^mH , (l) 

the core diameter of a vortex at position a? is on the 

_^ 1/2 

order of the local healing length £ — [gp (x )] ' , where 
AI is the particle mass, g is the interaction strength, and 
p = |\I/| 2 is the condensate number density, which varies 
smoothly over the sample from zero to its maximum 
value. This maximum is fixed by the chemical poten- 
tial p, which is not a parameter in the Gross-Pitaevskii 
equation, but an integration constant: the general solu- 
tion admits a time-dependent prefactor in ^ of the form 

Throughout this paper we will consider a quasi-two- 
dimensional condensate, so that g is dimensionless [ p6[ . 
Our results may also be applied, with trivial changes, to 
the case of a hydrodynamically 3D condensate in which 
the density is constant along the third direction, and 
the vortex lines are parallel to it. In three-dimensional 
background clouds, however, vortex lines generally bend 
[H and this introduces complications which usu- 

ally prohibit accurate analytical treatment J8| ; hence our 
restriction, in this paper, to the two-dimensional case. 
Even in the limit of a quasi-two-dimensional conden- 
sate, however, the inhomogeneous background density of 
a trapped condensate significantly affects the behaviour 
of two-dimensional 'point vortices'. The Thomas-Fermi 
(TF) surface, which is the edge of the condensate cloud in 
the hydrodynamic approximation p7[ , also affects vortex 
motion. 

It is typically the case, in current experiments, that 
near the vortex the trapping potential V varies slowly 
on the healing length scale. This small ratio of scales 
has two distinct implications for vortex physics. Firstly, 
it means that within many healing lengths around the 
vortex center, the gradient in V can be treated as a per- 
turbation. Secondly, it means that beyond a few healing 
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lengths from the vortex center, we may solve (|j) in the 
hydrodynamic approximation. These are two very dif- 
ferent kinds of approximation, but neither is valid every- 
where, and so they must be combined. Fortunately there 
is a significant intermediate region near the vortex within 
which both approximations are valid, and this will allow 
us to match the solutions they yield smoothly together, 
to obtain a complete solution. This is an example of a 
general procedure, known as boundary layer theory, or 
the method of matched asymptotics. 

Boundary layer theory is more than just a technical 
trick without physical meaning: it is an essential feature 
of their physics that vortices are small, healing-length- 
scale structures, realized within a larger-scale hydrody- 
namic medium. Boundary layer theory is a very direct 
expression of the multiple-scale nature of vortices, and 
the only alternatives to it are exact (or numerical) so- 
lutions, which will of course also exhibit the problem's 
multiple scales, or the replacement of some part of the 
boundary layer analysis with some uncontrolled approx- 
imation. The boundary layer method has yielded gen- 
eral equations of motion for dark solitons in quasi-one- 
dimensional inhomogeneous backgrounds p8|| , and has 
already been applied to 2D vortex motion in inhomoge- 
neous backgrounds 0, In this paper we will take 
a precisely similar approach. Our main contribution 
is an exact solution to the outer, hydrodynamic part of 
the problem, in a physically relevant and nontrivial case. 
Comparing our exact solution to the assumptions of pre- 
vious work will illustrate some general aspects of vortex 
physics that may not have been fully appreciated hereto- 
fore. 



Organization 

This paper is organized as follows. In the following 
Section II we review the hydrodynamic approximation, 
and then identify a loophole in previous derivations of a 
general solution supposed to be valid in the neighbour- 
hood of a vortex. We show that in fact a nontrivial 
problem must be solved, which generically lacks small 
parameters and involves the whole condensate globally. 
We then present a physically relevant case in which this 
problem can be solved exactly, in terms of a special func- 
tion which is only moderately obscure, and whose asymp- 
totic behaviours close to and far from the vortex can be 
obtained analytically. 

In Section III we follow earlier authors in obtain- 
ing the inner solution, perturbing around the numerical 
solution to the Gross-Pitaevskii equation for a vortex in 
constant background density. We then go on to match 
the inner and outer solutions together, and thereby de- 
termine the velocity at which the vortex moves parallel 
to the TF surface. With these results, we compute the 
free energy of a vortex in a moving frame, and thus assess 



the velocity at which the energetic barrier preventing vor- 
tices from entering the condensate will disappear. We 
compare these results to recent analogous calculations 
for critical rotation frequencies of harmonically trapped 
condensates, concluding that the latter overestimate crit- 
ical frequencies of large condensates by factors of order 
unity. 

In our final Section IV we discuss our results, interpret- 
ing the vortex motion as due to vortex buoancy, through 
a Magnus effect which is renormalized by the distortion of 
the flow field. We interpret this distortion as an 'infrared 
dressing' of the vortex, and emphasize that such dressing 
is a much more general effect than can be described with 
image vortices. We argue that the energetic barrier to 
vortices actually disappears at or below the surface mode 
critical velocity |23|, so that, strictly speaking, vor- 
tices entering condensates should not be said to 'nucleate' 
(cross an energy barrier by thermal fluctuations or quan- 
tum tunneling) . We then conclude with a brief summary 
and outlook. 



OUTER SOLUTION: HYDRODYNAMICS 
Hydrodynamic approximation 

The hydrodynamic approximation works as follows. If 
we define '3/ = ^fpe %B for real p, 9, the Gross-Pitaevskii 
equation becomes exactly 



dtp = ■ pS/9 (continuity) 



dt0 = - 



h 



Po 



M 
P{x ) 



M, 



V - 



e Po 



(2) 
(3) 



where the vortex center is located at xq. We write g in 
the form 1 / (£ 2 p ) in order to introduce the vortex scale £ 
explicitly. We seek a solution in which the only time de- 
pendence is due to the vortex motion, so that sufficiently 
far from xq we can set p — 0, and 9 = —p/H with con- 
stant p. If 9 and V are functions of ex/t; for some small 
e, then the second equation (||) just yields the Thomas- 
Fermi density profile Ptf(%) = M (p — V) / (ft 2 .g), up to 
corrections of order e or higher. The surface on which 
p TF vanishes is the TF surface, where the TF approxima- 
tion breaks down. (An additional boundary layer treat- 
ment is therefore required very close to the TF surface 

In this case we may expect (and later confirm) that the 
vortex velocity v vtx is order ec, where c = H/ (M£) is the 
speed of bulk sound near the vortex. This means that 
the corrections are formally of order e 2 . But there will 
also be corrections that diverge as one approaches the 
vortex, the leading one being the vortex kinetic energy 
cx r~ 2 . To keep these corrections no larger than order e, 
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we can only use the hydrodynamic approximation in an 
'outer zone' whose inner boundary is a circle around the 
vortex centre, having radius R of order e -1 / 2 £. In the 
'inner zone' inside this radius, we will use perturbation 
theory on the potential gradient instead. As long as it 
has the right order of magnitude, the precise value of R 
is arbitrary: it is merely a bookkeeping device to let us 
merge two different approximations, and all i?-dependent 
terms necessarily cancel when the two zones are patched 
together. So the full calculation will in fact be accurate 
up to corrections of order e 2 after all. And this definition 
of the outer zone also formalizes the requirement of being 
'sufficiently far' from the vortex to set p = and 9 = 
—fi/h, since for \x — xq\ > R the corrections v v tx 1 Vp, 9 
will be of order e 2 . 

The higher order corrections can all be computed triv- 
ially once the phase field 9 is known to zeroth order. 
(We will not actually compute such corrections in this 
paper.) The zeroth order phase field can be obtained by 
solving the continuity equation (^) with the zeroth order 
TF density, and so this is the main problem of the outer 
zone. 



Dual fields 

Because the phase will be multi-valued in the presence 
of vortices, it is convenient in two dimensions to define 
the field A which is dual to the velocity field (in the sense 
in which electric and magnetic fields may be dual): 



(A x ,A y ) = (d y 9, -d x 9) 



(4) 



The continuity equation then becomes V x (pA) = 0, 
which can be identically satisfied by setting 



A= 



(5) 



for some potential F(x,y). We will refer to F as the 
dual potential to the phase 9, but it is also known as 
the streamline function, because its contours of constant 
height are streamlines of the velocity field. Because it 
is single- valued even when vortices are present, F is a 
more convenient substitute for 9. F is not arbitrary, 
because we can see from (^) that V ■ A = V x V9 = 
2ttS 2 (x — xo) (since a singly quantized vortex is located 

atx ). Hence F must satisfy p V- (^p~ 1 VF^J = 2nS 2 (x— 

xo). But since we are only seeking the zeroth order phase 
field, we can replace p — > p TF in both (|^) and in the 
constraint on F, writing A = {po/ Ptf) as we H as 



VF 

Ptf 



2ttS 2 (x — xo) 



(0) 



where the = means the neglect of terms that will ulti- 
mately lead to corrections of order e 2 . (Introducing o 



subscripts on F and A would too greatly encumber our 
notation.) Eqn. (^J) is the essential equation determining 
the zeroth order outer zone phase field of a vortex. Note 
that since it is a linear equation, solutions with more 
vortices can be obtained trivially from the single vortex 
case. 

As a simple example, consider the classic case of a 
vortex in a constant background density profile, near a 
hard wall. Choose the y axis to run along the wall, and 
let the vortex be located on the x axis at the point xo- 
Since in this case p/p = 1, we merely have a Laplace 
equation to solve, with a delta-function source, and the 
boundary condition that there be no flow through the 
line a; = 0. Since adding a constant to F obviously 
does nothing, we can enforce this boundary condition by 
setting F(0, y) = 0. This is satisfied by the solution 



xo? + V 2 



1 (x 

Fhw — 7; m 7 — ; vT" 1 — 2 

2 (x + xo) 1 + y 2 



where a possible extra term proportional to x must van- 
ish in the case that the velocity field vanishes at infinity. 
This solution can evidently be obtained by the method 
of images, with the denominator of the logarithm repre- 
senting the dual potential of an image antivortex located 
at — xq- 

Obtaining an analogous solution for an inhomogeneous 
p is generally much more difficult, however. 



general results 



General equations of motion for point vortices in in- 
homogeneous superfluid backgrounds have nevertheless 
been presented, having been derived using the bound- 
ary layer approach [Q, [TTJ . According to these proposed 
equations, vortex motion is determined by the local trap- 
ping potential in the immediate neighbourhood of the 
vortex (that is, can be given by a universal formula in- 
volving the potential, its gradient, and its Laplacian at 
xo)- We pause here to point out a loophole in the deriva- 
tion of these results, which unfortunately allows correc- 
tions to vortex motion that are comparable in size to 
the proposed general results, and that in general cannot 
be computed without solving the global hydrodynamic 
problem. 

The approach of Rubinstein and Pismen, and of 
Svidzinsky and Fetter following them, is the same ap- 
proach we are following, of matching a hydrodynamic 
outer zone with a perturbative inner zone. At the point 
we have now reached, namely solving Eqn. (^|), these au- 
thors employ the variable H(x,y) = y p /pF(x, y) (Ref. 
p"T[ uses a more elaborate notation for the same function) 
and so obtain the equivalent equation, 



[V 2 -k 2 ]H = 2%8{x- x )S(y) 



(7) 
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for fc 



V p 



Vp 



-K n 



Their general procedure, for any slowly varying po- 
tential, is then to replace k(x) — > k(xo) = fco, 
so that (0) becomes simply a 2D Helmholtz equa- 
tion (with a delta- function source). It is then ob- 
served that the modified Bessel function H(x,y) — 

kQ\J{x — xo) 2 + (y — j/o) 2 ^ is a particular solu- 
tion to (0) with k{x) — > fco- Of course only the 
(x, y) — > (xo,yo) limit of this result is to be taken se- 
riously, so the Rubinstein-Pismen result for the outer so- 
lution near the vortex is 

Frp (x + £r cos cj> , y + £r sin (j>) 

= ^l+^r-Vp, ( 7 + In + 0(e 2 ). 

Euler's constant 7 = 0.577 appears through the asymp- 
totic behaviour of the modified Bessel function at small 
argument. 

The problem with this computation is that it assumes 
that the solution for H contains no components regular 
at (x,y) — > (xa,ya) (terms like 7q , 7i cos </>, etc.). Ac- 
tually, nothing ensures this. And so all that can really 
be claimed, on the basis of purely local analysis near the 
vortex, is that 

( i~ A \ 
lim F = 1 H — r ■ Vp I lnr 

X^XQ \ 2 J 

+k £r(B y cos (j) - B x sin 0) + 0(e 2 ) (8) 

for some undetermined constants B XiV . (A constant term 
can be dropped because it will not contribute to the dual 
velocity field A.) This unknown B implies an unknown 
correction to the superfluid velocity field near the vortex, 
and hence to the vortex velocity. So the strictly local 
approach, based on solving (Q) with k(x) — ► fco, does 
not actually lead to any conclusion at all for the vortex 
motion. 

In Reference M , Rubinstein and Pismen are careful to 
state that they are deriving the vortex velocity relative 
to the 'ambient' superfluid velocity, thus recognizing the 
corrections they are omitting. It is worth emphasizing, 
however, that the fluid velocity that may be called 'am- 
bient' is in a sense scale dependent. The effect omitted 
in Refs. j?], is indeed a flow near the vortex that is 
constant on the healing length scale; but it is generally 
not constant on the much longer hydrodynamic scale of 
the outer zone. So while such a flow may be called 'am- 
bient' from the point of view of the vortex core, in general 
it is not simply the flow that is specified at infinity (or 
wherever boundary conditions are imposed). 

Whereas Rubinstein and Pismen recognize but omit 
them, Svidzinsky and Fetter argue in Reference Jll|] that 



additional terms regular at the vortex may be ruled out, 
because the /„ functions all diverge at large argument, 
and so would violate any physical boundary conditions. 
While it is true that the /„ all diverge at large argument, 
they only do this on the scale fco, and once kty\x — xq\ 
is no longer small, the approximation k(x) — + fco is no 
longer good. Hence the argument of [O] is invalid. 



Motivation for an exact solution 

The moral of this general discussion should be clear. 
Implementing the hydrodynamic approximation has ex- 
hausted the benefits of the slow variation of the trapping 
potential on the healing length scale and thus, in general, 
Eqn. (JsJ) contains no small parameters. So except in spe- 
cial cases where additional small parameters may appear, 
such as for a trap with a high aspect ratio, or a vortex 
very close to the center of symmetric trap, no accurate 
equation of motion for a vortex may be obtained without 
solving a nontrivial hydrodynamic problem exactly. 

In JtJ, Rubinstein and Pismen provide one nontrivial 
exact solution, for the case described in our variables as 
M\jj,-V{S)] = h 2 gp [I (r/K)}~ 2 for some K > f . They 
focus on this case because it ensures that k(x) in (|?j) re- 
ally is constant, so that their solution quoted above be- 
comes exact. The trap required to realize their solvable 
example with trapped dilute condensates is not implausi- 
ble (it closely resembles a Gaussian well) ; but one would 
also need a finely tuned chemical potential, to make the 
condensate density just vanish at infinity. And it is 
a special feature of this finely-tuned case that there is 
no Thomas-Fermi radius. (Corrections to the Thomas- 
Fermi profile begin to exceed (£/7?.) 2 for r > 7?. In (7£/£), 
but their onset is very gradual.) It is therefore the main 
contribution of this paper to provide, in the next sub- 
section, another instructive exact solution to the hydro- 
dynamic problem of a superfluid vortex in an inhomo- 
geneous background, which is more directly relevant to 
currently typical experiments. 



Vortex hydrodynamics in a linear density profile 

The plane linear approximation 

Near a Thomas-Fermi surface, the potential is approxi- 
mately linear; and the TF surface of a large condensate is 
approximately flat. This motivates considering the ide- 
alized problem of a linear ramp potential 23 , 24|, ^?], . 
Choosing the y axis to run along the straight TF surface, 
we have the TF profile 



P = Pq- 



Xq 



(9) 
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Since it will turn out that F decays on the distance scale 
of xq , the linear ramp potential will indeed be reasonably 
accurate for real traps, which are obviously not globally 
linear, as long as the vortices are not too far from the TF 
surface; for harmonic traps, this means that xq should be 
much less than the TF radius. On the other hand, with 
V = ji — Xx (obtained automatically by taking x = 
on the TF surface), we have £ = H/^/MXxq, and hence 
e = £/ x = fi/ \/ AIXxq. The hydrodynamic approxima- 
tion breaks down within a few surface depths of the TF 
surface, where the surface depth p7f is 



5 = (h 2 /2MX) 



1/3 



(10) 



We therefore also have e = \/2 (5/xo) 3 ^ 2 , and so our 
boundary layer treatment of the vortex will require xq 3> 
5. See Figure 1 for a sketch of our model system indi- 
cating the relationships between the various length scales 
that will appear in our calculations. 

Vortices located well outside the TF surface (xq <C —S) 
can be described perturbatively as surface excitations of 
the condensate p4| : but vortices with centers within the 
surface layer |a;| < 5 require a nonperturbative treatment, 
which we are unable to provide. In our final Section, 
though, we will argue that nothing remarkable can hap- 
pen in this regime, which is analytically intractable but 
physically trivial. As mentioned above, the TF surface 
requires a boundary layer treatment of its own , inde- 
pendent of the vortex; but the post-hydrodynamic effects 
within this layer can easily be shown to produce correc- 
tions smaller than all results we will report by factors of 
at least (<5/:eo) 2 . Hence for vortices many surface depths 
into the condensate, we will be able to ignore completely 
the breakdown of hydrodynamics at the TF surface it- 
self. 



The equation 

Setting yo = by our choice of origin, Eqn. (^|) in this 
case becomes 

[dxx - ^<9x + d yy ]F = 2nS(x - x )S(y). (11) 

As warned above, this equation has no small parame- 
ters in it: even xo may be scaled away. Fortunately, 
however, (11) is exactly solvable. Taking advantage of 
translational symmetry in y, we can write 



/oo 
dkf k (x) i 
-oo 



,iky 



to obtain the ODE 



fk ~ ~fk ~ k 2 f k = S(x - x Q ) 



(12) 



(13) 




FIG. 1: Sketch with scales. Darkness of the background in- 
dicates condensate density, increasing linearly with distance 
x from the Thomas- Fermi surface at left, on which the y axis 
is placed. The hydrodynamic approximation is used in the 
Outer Zone, further than R from the vortex center; in the In- 
ner Zone the full Gross-Pitaevskii equation is solved with the 
potential gradient treated as a perturbation. The relative 
sizes of the four scales £,5,R and xo are indicated qualita- 
tively, but the ratios should obviously be rather greater for 
the analysis in the text to be accurate. The counter-clockwise 
circulating vortex is shown in the text to move in the negative 
y direction. 



We begin by solving the homogeneous equation, setting 
the RHS to zero. Multiplying f k by x allows us to rec- 
ognize a modified Bessel equation, and so we obtain so- 
lutions involving modified Bessel functions of order 1: 



f k = x[A k K x {kx) + B k I x {kx)] 



(14) 



for constants A k ,B k - We wish to impose the hydrody- 
namic boundary condition of no velocity through the TF 
surface, corresponding to F(0, y) = a constant that we 
can choose to be zero; so since \\m x ^Q K\{x) — l/x, we 
need A% = 0. We can accept no solutions that grow ex- 
ponentially as x —>■ oo, and so we must also set all Bk = 0. 
The only exception that might occur is the special case 
k = 0, in which we might under some circumstances re- 
quire /q = Bqx 2 , since F oc x 2 corresponds merely to a 
constant velocity field in the y direction. But for the case 
where the velocity vanishes at infinity, we need Bq = 
as well. 

Note the difference here between our rejection of addi- 
tional homogeneous solutions, and their rejection in [pj: 
we have solved our hydrodynamic problem exactly for 
all x, and not just asymptotically near the vortex. Of 
course, to apply our model to a realistic case we must ad- 
mit that the potential becomes nonlinear at sufficiently 
large x or y; but as long as this nonlinearity scale (call 
it A) is much larger than xq, our exact solution will be 
applicable over an essentially infinite region, as far as the 
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vortex is concerned. One loophole does remain, however, 
even in this case. We have fixed Bq = by demand- 
ing that the velocity field vanish at infinity, to which we 
assume that our linear model of the potential extends. 
But in a real, finite system, even if it is modelled well 
by a linear potential over a large region, it could be that 
boundary conditions farther away from the vortex than 
A imply some nontrivial, nonzero Bq. This would allow 
corrections to the vortex velocity that one can expect to 
be of order h/(AM). To compute them would require 
solving the actual, finite hydrodynamic problem. 

The vortex solution 

Having solved our homogeneous equation in general 
and concluded that the only solution admitted by our 
boundary conditions is zero, there remains the main task 
of obtaining a particular solution for F with the delta 
function source at x = xq. Since we have the general so- 
lutions to the homogeneous equation, we can patch them 
together to meet our boundary conditions and also fit the 
delta function source, by writing 

_ x_ ( K 1 (kx )I 1 (kx) , x < x 1 
/fe ~ x \ hikx^K^kx) , x > x Q J 

for a constant C. We then fix C by imposing that the 
discontinuity in the derivative of fk at x = x be equal to 
1, as required by @. Since KKOMO " ^WiCO = 
— 1/£, this gives C = — 1, and so we have 
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FIG. 2: Hydrodynamic flow lines around a 2D point vortex 
near (left) a TF surface and (right) a hard wall. In both 
cases the surface of the condensate is the vertical axis, with 
the axes marked in units of xq. The left plot is constant 
height contours of F = — \J x j x§Q\i-2,(z) for z(x,y) defined in 
the text. The right plot consists of constant height contours 
of the homogeneous dual potential with a hard wall surface, 
Fhw = MVfc - xo) 2 + y 2 /y/(x + x ) 2 + y%. 



The vortex flow field can be visualized by plotting 
contour lines of constant F which correspond to fluid 
flow lines. Such contour plots are available numerically 
through common commercial software applications that 
include large libraries of special functions. An example is 
shown in Figure 1, with the corresponding flow pattern 
for a vortex in homogeneous density near a hard wall, 
for comparison, ft is obvious that the density gradient 
distorts the flow pattern signficantly. 



F = —2 — / dk cosky lKi(kxo)Ii(kx)9(xQ — x) 
xo Jo 

+K 1 (kx)h(kx Q )9(x - x Q )} (16) 

where 8(x) is the step function. Happily, this integral 
can be evaluated in closed form in terms of a Legendre 
function of the second kind [M, yielding 



where z 



F(x,y) = -y/x/x Qi (z) , 
— , — ) is defined as 

(x - xq) 2 + y 2 



y 



2xxq 



= 1 



2xxq 



(17) 



(18) 



The Legendre functions of order v satisfy the Legendre 
differential equation 



d 

dz 



1-z 2 



dQ v 



dz 



= -v{v+l)Q v 



(19) 



and the functions of the second kind have logarithmic 
singularities at z — ±1. It is straightforward to ver- 
ify, using (|l9|), that ( p7[ ) satisfies (^]) with the boundary 
condition F(0, y) = 0. 



Outer asymptotics 

At large z, Q 1/2 -> 2- 5 / 2 nz~ 3 / 2 (easily obtained from 
(^B|) below, or see [pi)), so that the dual potential F falls 
off at long range much more quickly than the long-range 
logarithmic behaviour found at constant background den- 
sity. Inserting our definition of z we therefore have 



lim F(x,y) 

\y\— >oo 



lim F(x,y) = -—— 



lim F(x,y) = 



TTX 2 Xq 1 


1 + - 




2x \ 




X 2 






2x1 V 


+ y - 


2 
Xq 



2\ -3/2 



2\ -3/2 



(20) 

(21) 
(22) 



The kinetic energy density due to the vortex thus falls 
off as x~ 5 at large x for fixed y, and as xy~ 6 for large 
y and fixed x, which allows us to conclude that vor- 
tices near Thomas-Fermi surfaces are essentially local- 
ized structures that do not extend into the condensate 
beyond a depth of order xq . 



7 



Inner asymptotics 

For matching with the inner zone solution, which will 
determine the vortex velocity, we need the behaviour of 
Qi/2( z ) a s z — > 1 + . This can be obtained analytically. 
Using dimensionless polar co-ordinates (r, <j>) centred on 
the vortex, 



x — x + cos <j) , y — £r sin <fi, 



(23) 



and taking e = £/%o, we can express z as we approach 
the vortex as 



z (1 + ercos 4>, ersin</>) = 1 + 
= 1 + 



e 2 r 2 
2 1 + er cos < 
e 2 A 2 



(24) 



We can then use this expansion of z in the integral rep- 
resentation p2j of the Legendre function 



ds 



o (2 + \/ z 2 — 1 cosh s) 



—1 (25) 



to obtain 
Qi(l + £ 2 A 2 /2) = 0(e 2 ) + 



ds 



(1 + eAcoshs) 



3/2 



(26) 



Then we can note that 



ds 



(1 + eAcoshs) 
2 ds 



3/2 



(1 + eA cosh s 



,3/2 



1 + e s+ln ■ 



-fl+ln - 



3/2 



(27) 



and evaluate each of the two integrals using different ex- 
pansions in e (thus evaluating ( pq ) using a kind of mina- 
ture boundary layer theory of its own!). That is, in the 
first integral of the RHS of (^) we simply expand in eA 
and integrate term by term; and in the second we use 

ds 



-s+ln ■ 



ds 



In ^5 (l + e») 



3/2 



In 



1 



(28) 



in ■ 



where in the second line, as hereafter, = means dropping 
a term of order e 2 . 



Combining both results we obtain 



)i(l + £ 2 A72) 



2 + In 



eA 
~~8 



, er sr cos i 
- 2 - ln ¥ + — 



(29) 



and hence 



F(xq + C r cos 4>i ^i" sin < 



1 



er cos < 



1 £r 
In — 



-2 + -r cos 4> + 0{e 2 ) 



(30) 



Translating from the dual potential F back to A ac- 
cording to (pp, to first order in e we therefore have 



lim A r 

X — >XQ 



1 e cos 6 , er 
- H In — 

r 2 8 



lim = -- 

z^£ c; 



|(l + lnf) Mn, 



(31) 
(32) 



Applying (0) we then see 



ercos0 er 
Jim d^fl = 1-1 In — 

X^Xq 2 8 

£ / er \ 

lim d r 8 = — I 1 + In — ) sin (j) 

x — >xq 2 \ 8 / 

from which we can obviously extract the condensate 
phase field near the vortex, 

0(r,0)=0 + er^ln^ + O(£ 2 ). (33) 

Note that there is no trace of Euler's constant in this 
equation. 



INNER SOLUTION: GROSS-PITAEVSKII 
The perturbative problem 

In the inner zone we work entirely in terms of the di- 
mensionless polar co-ordinates centered on the vortex: 
x = xq + £rcos</>, y = £rsin0. (Since the actual con- 
densate density at the vortex center vanishes, it is per- 
haps worth clarifying that for the local healing length 'at' 
the vortex £ = (gp Q ) we use for p the background 
Thomas-Fermi density extrapolated to the vortex loca- 
tion.) To obtain our inner region solution, we expand 
the potential to first order about the vortex position: 
V(x) = Xx = \x Q [l + er cos + 0(£ 2 )] . It will then be 
convenient to rescale the wave function, defining 



efict 



(34) 



8 



where c = H/ (M£) is the local speed of sound at the 
vortex, and v vtx — eftc is the vortex velocity, as yet un- 
known. (We have reduced the number of symbols to be 
defined by anticipating the fact that the vortex velocity 
will be order e.) Then writing ip = ip + stpi + ... , we 
find 



1 



1 



c) 2 



i> + [Kl 2 - 1] % = (35) 



r 



f 



1]/ 



(36) 



as our zeroth order equation. The only solution which 
corresponds to a singly quantized vortex at r = is ipQ = 
/(r)e 4 * where /(r) can be taken as real, and satisfies 

1 
2 

and f(r — > oo) = 1. Note that a vortex velocity of 
order e° would require a v v t x ■ VV>o term on the RHS 
of (H|), and the only well-behaved vortex solution with 
such a term present is e lMv ' x / h if> ; but this would imply 
a phase gradient of order e° as one approaches the outer 
zone (r — > R/S,), and this would be inconsistent with our 
outer solution (|33"|). 

At first order, though, we find 

1 

~2 



1>i 



r cos <f)ip 

+ {0y COS 



COS ( 



13 y sin 



[2/ 2 - 1] 



r 



(37) 



It is not easy to solve this equation for tp^} but we are 
actually only interested in two pieces of information. We 
need to know the asymptotic behaviour of tp 1 many heal- 
ing lengths away from the vortex center, so that we can 
smoothly match the inner solution to the outer. And we 
need to determine /3, which will be fixed by the require- 
ment that does not blow up as r increases. Following 
Rubinstein and Pismen , we will be able to obtain this 
information without solving (|3^) explicitly. 
At large r we can write ip 1 — ► e 



^[S + iT] and find 
at large r, the 



S — ► | cos 4>. Then since / — ► 1 — (2r) 
leading terms in T are driven by the —ir~ 2 d c / ) S crossterm, 
giving the asymptotic equation V 2 T = r _1 sin0, with 
the solution 

T — ► r (a x cos <f> + a y sin ( 



— In r sin < 
2 



(38) 



for coefficients a that will be fixed by matching with the 
outer zone. (We drop a constant which can obviously be 
absorbed in tp , and which is of no consequence anyway. ) 
So we have 

T 

- (cos (p ■ 



lim * 



i In r sin 0) 

+iea ■ f] 
+ (er/2) (lnr + 2a y ) sine 
-era x cos d>] 



Vpexpi U 



r er . 
e-ln^sm, 



(39) 



dropping terms of 0(e 2 ). In the last line we have im- 
posed matching with the outer solution (B3f) to fix 



a x = 



v„ = — m — . 

y 2 8 



(40) 



To constrain p, note that differentiating (|35| ) with re- 
spect to x or y shows that d x ip and d y i/j are two inde- 
pendent solutions to the homogeneous equation for 
Writing £ as an abbreviation for the left-hand side of 
(|37|), and J as an abbreviation for the right-hand side, 
we integrate both sides of ( |37| ) with d x ipQ out to the large 
dimensionlcss radius Integrating by parts and us- 

ing our results for tp 1 at large radius reveals 



Re 



rdr <z> d<fi 



, "' [ f cos (p + i- sin cp ) t 



-Imfd^ [e~^ sin<t> (d r ^ + r~ Vi)] r= 



r=R/£ 



TV 


, £-R 1" 


7T 










2 


ln 8f + 2. 


2 










Re 






e 


















'J 


dr//' [r 


2 + 2/3 




D 







f cos </> + i— sin ) J 
r 



<ir r/ 2 



(41) 



This allows us to obtain 

'l+ln(E/0 



a 



lim 

R/£-*oo 



1 + 111 



i? 



ln/'(0) 
0.114 . 



dr 



1 







/ dr r(l - 


■-/() 






r 




f 


/"■ 

r — 





(42) 



In the second line we have used (|36|) to re- write (1 — 
/ 2 ) in the integrand, and in the last line the numerical 
evaluation comes from a numerical solution for /(r) . 



Similarly, using the dyijj 
solution, we find 



/3 X = 0. 



(43) 



(The numerical result in p2| ) was obtained with Math- 
ematica ]33] |, using two different methods of numerical 
solution for /(r) (shooting and relaxation), whose re- 
sults agreed with each other. Quite stringent settings 
of the options for starting step sizes, etc., were required 
to obtain this agreement, and in both methods the sin- 
gularities at r — had to be regulated, for example by 
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replacing r — > \Jr 2 + 10~ 20 in the singular co-efficients in 
( |36| ) . Our result does not quite agree with the evaluation 
reported in 0: their value of 0.405 for their quantity 
lnai corresponds to a value of —0.126 for our quantity 
(3 y - a y .) 

Combining ( f4(i| ) and (0), we obtain the total vortex 
velocity. It is in the negative y direction: parallel to the 
TF surface, and in the direction of the fluid flow between 
the vortex and the surface. Extracting explicitly the xo 

dependence hidden in e = \J 2S 3 /xq, we can express the 
magnitude v v t x of the vortex velocity purely in terms of 
xo and surface parameters: 



l - 
0.8 
0.6 
0.4 
0.2 



10 



12 



14 



Vvtx 



sc 



where 



J__!l 



1/5 1 



h 




(2,2x1 


j +0.114 


Mxo 


V * 3 , 


5 

V c 

xo 


>( 




0.980 





V c = 



h 

~M6 



(44) 



(45) 



FIG. 3: Velocity of a vortex parallel to a Fhomas-Fermi sur- 
face in a quasi- two-dimensional condensate (upper curve). 
Fhe vertical axis is velocity in units of the surface mode criti- 
cal velocity h/ (MS); the horizontal axis is distance of the vor- 
tex center from the TF surface, in units of the surface depth 8. 
As a leading order result in 5/xo, the curve is not meaningful 
to the left xo = 5. The dashed curves, presented for compar- 
ison, are the Rubinstein-Pismen result (short dashes) , and the 
velocity h/ (2Mxo) of a vortex the same distance from a hard 
wall in a condensate of constant density (long dashes). The 
sum of the two dashed curves is quite a good approximation 
to the solid curve. 



is the surface characteristic (and critical j24j) velocity. 
Eqn. J44| ) is the first main result of this paper; it is 
illustrated graphically by Figure 3. It will be accurate 
for vortices in quasi-two-dimensional condensates as long 
as the vortex distance from the TF surface Xq is much 
larger than the surface depth 8 but much smaller than the 
TF radius. For quasi-2D condensates of size comparable 
to current three-dimensional condensates, this will be a 
significant regime of validity. Of course, even apart from 
its realism, (Q) remains instructive as an accurate result 
for an idealized problem. 



In this expression for E we have subtracted off the 
Thomas-Fermi free energy of the vortex-free condensate. 
We will now evaluate this free energy to leading order in 
e using our results from above. As above, we will com- 
pute the inner and outer zone contributions separately, 
and add them. In this sum, all dependences on the inner 
zone size R ~ e _1 / 2 £ necessarily cancel out. 



Computing E 



Free energy and vortex penetration 

If dissipation occurs in a frame moving with respect 
to the condensate with velocity Vdis , then the vortex will 
drift towards or away from the surface, in order to mini- 
mize the free energy 



G 



TTh 2 



2Mg8 d 



7-n Vdis 

E — p 

Vc. 



(46) 



where the prefactor is a surface energy scale, so that the 
dimensionless energy E and y component of momentum 
p are given by 



Since we know that v vtx is of order e, it is not hard to 
see that the order e terms in the inner component Ei n of 
E vanish, and up to corrections of order e 2 we have 



/' 2 + ^ + (/ 2 -l) 2 



S 



dr d r 



rff' + r 2 f 



r2/ 4 + /2 _ r2/ , 21 



" S 



drr(l- 2f) 



xo 
6 



(49) 



E 



XQ 
TTp 5 2 



(I 2 :. 



Im / dxdy ^*d y ^. 



Fx 
9 



(47) 
(48) 



where in the last line we also drop terms of order 

which are not only of order e, but will also be cancelled 

by terms from the outer component E out . 

In the outer zone we must integrate the energy density 
over the entire half plane, except for the inner zone circle 
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of radius R ~ e 1 I 2 xq. Denoting integration over this 
region with the subscript A we have, to leading order in 

e 



E ut — —7 / dxdy — 
no J p 

A 



r 2 



71 6 



dxdy V 



V0 
F 



Xq 



:)•() 



= — / dxdy — 



.4 

XO 

71" (5 



-, F - 
T dS- -VF 

J TTO J X 
OA 



Fd r F 



r=R/£ 



(50) 



To obtain the second line from the first we use (|llj) , where 
the RHS vanishes everywhere in A. Within the second 
line we use the divergence theorem, where the surface 
terms at x = and at infinity all vanish, so that the only 
contribution is on the boundary with the inner zone at 
the circle of radius R about xq . So putting Ei n and E out 
together we have 



E 



_Xq_ 

S 

X 



( ln 



6 



0.980 



(51) 



Writing E as we have in the last line emphasizes that 
dE/dxo is proportional to v vtx , indicating that xq and 
v v t x are canonically conjugate, as one expects for a vor- 
tex. 



Computing p 

To compute p it is very helpful to note that the dual 
potential F can be extended into the inner zone, by using 
the continuity equation and the fact that the density is 
constant in the frame co-moving with the vortex to obtain 



= (M/h) [v vtx d y p - p] 

-> \ ( - M 
= V • p V0 + -ryv vtx 



which allows us to define 



pdyB = p d x F + —v vtx (p - —Po 



pd x 



-p d y F 



(52) 

(53) 
(54) 



Taking this definition of F into the outer zone, it clearly 
agrees with our previous one up to post-hydrodynamic 
corrections. (Alternatively one could repeat the outer 
zone analysis using this slightly different definition of F, 
and confirm that no differences arose up to order e 2 .) 
Furthermore one can show from the Gross-Pitaevskii 



equation that this extended F is perfectly regular as 
r — > 0, where it behaves as J drr~ l f 2 . We can therefore 
write 

p = — % / dxdy —d y 9 
ttS J Po 

dxdy (—-—) (55) 
\Po x oJ 

F(0,y)- lim [F(x,y)\ 



xo v vtx 



— \ i dy 

TTS J-00 



The leading contribution to p is the last line of 
evaluation is very simply obtained from (|2l|): 



Its 



X 



lim 



7tS x ^°°J- 



26 2 



2S 2 



lim 



dyF(x,y) 



d vf 1 + y_ 

X \ X 2 



-3/2 



dy 



00 (1 + ?/) 



2^/2 



(56) 



This leaves out the density deficit integral in the second 
line of (|55|), which can be shown to be of order S/xq, and 
hence smaller than the leading term by a factor of e 2 . 



Free energy curve 



We therefore have, up to corrections of order e 2 



G = 



irk 2 



2MgS z 



\ 2 I v c \ / 



(57) 

which is plotted in Figure 4. At least within the hydro- 
dynamic approximation, the energetic barrier to vortex 
penetration has clearly disappeared for \vdis\ > v c . This 
is the second main result of this paper. Since we are con- 
sidering a vortex with counter-clockwise circulation, lo- 
cated along the positive x axis away from a TF surface on 
the y axis, it is to be expected that a negative Vdis makes 
vortex penetration become energetically favourable. 



Application to rotating harmonic traps 

For a harmonically trapped condensate of spatial size 
Rtf, the nonlinearity of the trapping potential near the 
TF surface, and the curvature of the TF surface, can both 
be neglected in the limit where xq/Rtf — * 0. Hence our 
results can be applied to harmonically trapped conden- 
sates if we regard them as leading order approximations 
in both 5/xq and xq/Rtf- If we consider a rotational 
symmetric trap, we take Rtf to be the Thomas- Fermi 
radius in the plane of symmetry, and then assume ei- 
ther a quasi-two-dimensional 'pancake' trap, or a quasi- 
cylindrical 'extreme cigar' trap in which vortex lines are 
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FIG. 4: Free energy G/G 
Tvh 2 /(2MgS 2 ), 5 = (2MX/H 2 



-- E 

-1/3 



- PVdi 

and v, 



,/v c , where Go = 
= h/(MS). Hor- 
izontal axis is Xo/5, distance of vortex center from Thomas- 
Fermi surface in units of surface depth 8. From uppermost to 
lowermost, the curves are for Vdis/vc = —0.4, —0.6, —0.8, —1. 
The plots are shown dashed for xo/5 < 3, and stop at 
xo/S — 1, to reflect that the curves are leading order results 
in 5/xo- 



parallel to the long axis z. In the latter case the z length 
of the condensate will only appear as an overall factor in 
the free energy, and so effectively these two cases are 
exactly similar. We can translate our parameters into 
those commonly used for harmonic trapped condensates 
by noting that the potential gradient at the TF surface 
is just 



A = MuP-Rtf 
where u> is the radial trap frequency. This yields 

. , ft 2 ^ 1/3 
o 



2M 2 oj 2 Rtf 



(58) 



(59) 



with calculations. The prospect that an initially straight 
vortex line might bend more and more as it moved seems 
difficult to rule out.) 

We can also obviously interpret Vdis = QRtf where 
51 is the frequency at which a perturbation to the po- 
tential is rotated in order to stir the condensate. With 
this translation of Vdis, we can compare Figure 4 to the 
right-hand edge of Figure 5 from Reference jl3|], and 
compare ((57)) to Eqn. (49) of Q ('SF49') in the limit 
where £ — > 1 — 2xq/R in that equation. An essentially 
similar equation is presented in |l4| as their Eqn. (9). 
(Both these works examine three-dimensionally harmonic 
traps, in which, however, vortex curvature is neglected, 
so that the length of the vortex line is proportional to 
x\j 2 . Hence in their cases energy and momentum scale 

as Xq ^ 2 and x^ 2 , respectively, rather than xq and Xq as 
in ours.) Equation SF49 was derived by assuming that, 
for a rotationally symmetric trap, the phase field of the 
vortex would be well approximated throughout the con- 
densate by a simple e 1 ^ with <f> the usual polar co-ordinate 
centered on the vortex. Unless the vortex is very close 
to the center of the trap, this ansatz violates continuity 
significantly over most of the condensate. And we know 
from our results above that it significantly exaggerates 
the degree to which the velocity field of a vortex cen- 
tered near the TF surface extends into the bulk of the 
condensate. Hence according to SF49, the rotation fre- 
quency Vly [Q at which the free energy becomes negative 
for all xq > 5 is an overestimate (i.e. is greater than 
v c /Rtf) by the significant factor (5/4) \u(Rtf/^o)- So 
while discussions of vortex free energies based on the sim- 
ple e 1 ^ ansatz will typically be qualitatively sound, they 
can easily be inaccurate by factors of two or more, and 
cannot be used to obtain precise predictions for critical 
velocities or stirring frequencies. 



and 



Mu^JxIRtf 



3/2^1/2 
x ll TF 



(60) 



for the trap size ao = (Mlo/K) 1 . If we estimate the 
size of e in an effectively 2D condensate by supposing di- 
mensions typical of current large three-dimensional con- 
densates, taking ao ~ (im and Rtf ~ 25/im gives e < 0.2 
for xo > ao- So we see that the calculation should be 
accurate for vortices a few microns inside the TF surface. 
(It should be emphasized once again that three dimen- 
sional effects such as bending of vortex lines are known to 
be present in current experiments, but are entirely absent 
in our model. It seems reasonable to hope that vortex 
bending will be slight and have small impact in highly 
prolate or highly oblate traps, which should approach ei- 
ther our 2D or cylindrical limits; but it is difficult to ex- 
tend the present analysis far enough to support this hope 



DISCUSSION 

Magnus effect and infrared dressing 

The vortex velocity component /3 — a y is perhaps the 
least trivial aspect of the vortex motion, inasmuch as 
it describes a component of the motion that is deter- 
mined by the inner zone analysis alone, independent of 
the background fluid velocity a y in its immediate neigh- 
bourhood. (The latter is nontrivial to determine, from 
the outer zone hydrodynamics, but trivial in the way it 
moves the vortex). As first pointed out by Rubinstein 
and Pismen [Q , this intrinsic motion is along contours of 
constant Thomas-Fcrmi density. We can understand it 
qualitatively by noting that the vortex is a bubble-like 
density defect, and so experiences a buoancy force in the 
opposite direction to the trap force. Due to the Magnus 
effect which dominates vortex motion in superfluids and 
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normal fluids alike, this force produces not an accelera- 
tion, but a velocity at right angles: 



F x z F x z 



VMag = 



Mnp 2nhp 



(61) 



since in our case the vortex circulation k is 2ttH/M. (If 
we had taken a vortex swirling in the opposite direc- 
tion, we would indeed have found (3 y to have the opposite 
sign.) 

The precise magnitude of the buoancy force, and hence 
of the vortex velocity, is nontrivial, because the naive 
buoancy force 

F = VVJ d 2 x (p - |*| 2 ) 

= -XEc(2nph) ( drr(l-f 2 ) (62) 



is logarithmically divergent. What fjj showed first, and 
our own inner zone analysis has reproduced, is that the 
buoancy force is renormalized. The longer range effects 
of the potential gradient produce a logarithmic distortion 
of the vortex flow pattern, and since the classic Magnus 
effect is for a strictly cylindrical flow, this distortion pro- 
duces the counterterms (1 + lni?/£)/2 in the first line of 

On the other hand, we can see from ( |39| ) that the 
flow distortion in implies a velocity field component 
oc ylnr, which is essentially indistinguishable in any 
range of r from a background velocity perpendicular to 
the potential gradient. So the distinction between the 
intrinsic (/3 — a y ) and ambient a y components of vor- 
tex motion is somewhat artificial, and it is more natural 
to consider the whole flow pattern as a 'dressed vortex'. 

Previous analytic studies of vortices in BECs have of- 
ten neglected this dressing effect, by assuming that the 
phase field e l8 ^ of a vortex always remains e 1 ^ in polar 
co-ordinates centred on the vortex. Close to the core 
of a vortex this is indeed a good approximation, as long 
as 'close' means that the background condensate density 
has not varied appreciably. If further away from the 
core the background density varies (and is not rotation- 
ally symmetric about the vortex) , then thia simple ansatz 
for the phase field obviously fails to satisfy the continuity 
condition. Hence at distances from the vortex core that 
are on the scale of the trapping potential's variation, the 
phase field will depart significantly from e 1 ^ . Further- 
more, since the corrections to e 1 ^ are caused by density 
variations on the potential scale, their own spatial scale 
will be of this same order. And this means that the 
corrections will extend into the vortex core region, in the 
form a 'local ambient' flow. This flow is constant on the 
healing length scale, but it exists just because the vortex 
is present in the inhomogeneous sample, and in this sense 
can be considered part of the vortex: it is a component 
of the vortex's infrared dressing. 



Against image vortices 

The simplest example of this infrared dressing phe- 
nomenon is already well known, as it occurs in the spe- 
cial case of inhomogeneity that is a 'hard wall' boundary 
on an otherwise homogeneous sample. Here the non- 
trivial requirement of continuity far from the vortex is 
simply the constraint that there be no flow through the 
hard wall surface. In this special case the hydrodynamic 
problem to be solved for the vortex phase field is simply 
the Laplace equation, with Dirichlet boundary conditions 
on the surface. For many shapes of surface, the method 
of images is a technical trick that produces the required 
solution. This convenient technique has perhaps had 
the unfortunate side-effect of obscuring the general phe- 
nomenon of vortex infrared dressing, by giving rise to an 
impression that 'image vortices' are the only significant 
effect of inhomogeneity. In fact the method of images 
is restricted to conveniently symmetrical equations and 
boundary conditions, and it requires an exact solution 
which is known in the Laplace case, but not in general. 
It is of no help in the case of a Thomas-Fermi surface, 
because the single vortex solution automatically satisfies 
the correct boundary condition. And it is not even ap- 
plicable for a hard wall surface, unless the density profile 
within the wall has a fortunate form. (For example, the 
solution obtained in this paper will allow image solutions 
for a hard wall surface perpendicular to the gradient of 
a linear potential, but not for walls at other angles.) 

Thinking about image vortices will typically allow a 
qualitative understanding of how vortices will move near 
a surface; but this will not be quantitatively reliable. 
For those strongly attached to the image vortex picture, 
the singularity of Q v {z) at z — — 1, corresponding to 
x = —xo, may reveal an image vortex even in the case 
solved above. But since image vortices are actually just 
a mathematical trick, for representing the physical effect 
of a surface, it would arguably be just as well not to rely 
on them as a conceptual tool beyond their regime of real 
applicability. Instead of picturing repulsive or attractive 
image vortices, it is not so hard just to think about the 
accelerated flow through the 'channel' between a vortex 
and a surface, and so obtain a qualitative understanding 
that is equally convenient and more genuine. 

On the other hand, using the Laplace image vortex 
velocity h/(2Mxo) to estimate the size of surface effects 
on vortices is obviously better than nothing, if the full 
hydrodynamic problem is intractable (as it may well be) . 
Combining this crude image vortex theory with the lo- 
cal theory of Rubinstein and Pismen H does seem to 
work surprisingly well for the particular case analysed 
in this paper; but there are no grounds for expecting it 
to be generally accurate. Using such uncontrolled ap- 
proximations to guide experimental design might be rea- 
sonable, given the many unknown factors present in the 
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early stages of an experiment. But even fairly substan- 
tial disagreements between experimental measurements, 
and theoretical predictions based on such zeroth-order 
estimates, would be no evidence of novel phenomena. 

Estimating the critical velocity for vortex 
'nucleation' 

Since the vortex flow field in a homogeneous sample 
is already long ranged, it is possible that the infrared 
dressing which occurs in an inhomogeneous background 
may distort the long ranged velocity field by actually 
reducing its extent. We have seen in this paper that 
for a vortex near a Thomas-Fermi surface this sort of 
screening effect docs indeed occur, so that the vortex is 
a well localized structure. As we discussed briefly at the 
end of the last Section, this effect can be a significant 
correction on estimates of the velocity of stirring needed 
to drive vortices into a condensate. 

Since surface Bogoliubov modes can be considered as 
the motion of vortices in the ultra-dilute tail of the con- 
densate density profile extending outside the TF surface, 
the present paper can be regarded as a complement to 
p4j , showing how vortices may behave soon after they 
have 'nucleated' (that is, entered the condensate cloud). 
As we have mentioned, our hydrodynamic and boundary 
layer approximations are only valid if the depth of the 
vortex inside the TF surface is much greater than the lo- 
cal surface depth, and also much smaller than the sample 
size. (Hence the applicability of the plane linear model 
for vortices inside the condensate is somewhat more re- 
stricted than for the Bogoliubov surface mode spectrum 
in , where it is only required that the surface length be 
much less than the sample size.) Because of its inapplica- 
bility at small xa , Figure 4 is not really adequate to derive 
a precise result for the critical Vdis above which vortices 
will tend to enter the condensate spontaneously; but one 
can deduce from it with good confidence that this critical 
value must be above v c /2. The lowest curve of Figure 
4 certainly shows that, above the surface mode critical 
velocity v c , no barrier to vortex penetration will remain 
within the hydrodynamic region. This implies that the 
limitation on vortex penetration is in the perturbative 
region, where the surface mode analysis of pij fixes the 
critical velocity at v c — h/(M8). And this means that 
vortices enter a condensate through an ordinary insta- 
bility, and not by a quantum or thermal barrier-crossing 
process, which is usually what is meant by the term 'nu- 
cleation'. But two warnings must be attached to this 
conclusion. 

The first is that neither the surface mode analysis of 
pij nor the boundary layer theory of the present paper 
are valid for vortices within a few 8 of the TF surface, and 
so one might in principle worry that a narrow energetic 
barrier might still exist, within this region, at velocities 



above v c . Resolving this concern analytically would be 
very difficult, but it can be dismissed with physical ar- 
guments. The TF surface is not a physical surface, not 
a skin or a wall; it is a place where the condensate den- 
sity almost vanishes. In the neighbourhood of the TF 
surface one can expect to see a transition between in- 
ner and outer regimes, but there is simply not enough 
of anything there for this interface to constitute a third 
regime of its own. Furthermore, vortices whose centers 
are within a surface depth of the surface have core sizes 
on the order of 5 as well, so that a vortex centered on 
the TF surface already extends into the hydrodynamic 
regime on one side, and the perturbative regime on the 
other. This makes it very implausible that a barrier 
arises at the TF surface, when the free energy is mono- 
tonic in the same direction on both sides of it. And, 
finally, experience with numerical integration of dissipa- 
tive 2D Gross-Pitaevskii equations has always shown that 
once vortices begin sinking towards the TF surface, they 
pass through it and enter the condensate without any dif- 
ficulty. Hence we conclude that at velocities above v c no 
barrier exists for vortices, at least in condensates large 
enough for the linear density profile to be an accurate 
local model. 

The second warning is more important, which is that 
three dimensional effects may perhaps hold vortex loops 
near the TF surface, simply because for a vortex line to 
sink more deeply into a condensate it must typically grow 
in length, adding the vortex 'string tension' to the free 
energy. If the vortex line bends and loops, this effect 
might be much greater than can be accounted for by the 
factors of x^ 2 allowed in recent 3D calculations assuming 
straight vortex lines. Just how strong such an effect 
might be is far beyond the scope of this paper; but there 
have been some indications in experiments that stirring a 
condensate at just above the critical frequency generates 
twisting vortex loops that remain near the edges of the 
cloud @. 

Summary and outlook 

This paper has presented an exactly solvable hydrody- 
namic problem, namely the case of a point vortex in a 
plane linear background density profile, which can sup- 
plement the hard wall as an example on which to base un- 
derstanding of general cases of vortices interacting with 
surfaces. This particular solution is also directly relevant 
to the real problem of a vortex near the Thomas-Fermi 
surface of a large condensate. It will even be an accurate 
approximation in a realizable regime, because the veloc- 
ity field in this case actually falls off much more quickly 
than in the homogeneous case, so that the longer-range 
effects of nonlinear potentials and curved surfaces may 
be neglected more generally than one might initially have 
expected. 
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It is unfortunate, however, that the two conditions of a 
density gradient and a surface are combined in this solv- 
able case, because it therefore does not afford us much 
intuition about how their effects may differ. It should 
be clear from the preceding discussion that the effects 
found in this case are not purely surface effects: if at 
some point the linear trapping potential levelled off and 
then diminished, the TF surface would be replaced by 
a mere local minimum in background density, but the 
density gradient alone would still have effects on vortex 
dynamics. Vortices with density inhomogeneity far from 
any surfaces should be analytically tractable in some sim- 
ple models, such as slightly varying periodic potentials. 
A few aptly chosen numerical examples might be almost 
as instructive, and require less effort. 
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